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The variability inherent in rhythmic physiological pat- 
terns may arise from both stochastic and deterministic 
sources. Though variability is disruptive when in excess, 
it may also be functional in evolved systems. Here we 
focus on the neural control of respiration which is critical 
for survival in many animals. The sources of respiratory 
variability, as well as its possible function, are unknown. 
A fundamental component of respiratory pattern genera- 
tion is the preBotzinger complex (preBotC) which is part 
of the ventral respiratory column located in the brain- 
stem. Here we use in vitro and in silico models of the 
preBotC to study the variability of respiratory rate. 

The in vitro model exhibits a bounded range of variabil- 
ity in which the upper and lower limits are functions of 
the respiratory rate. The in silico model consists of two 
types of neurons: intrinsically bursting "pacemaker" cells 
and tonically spiking cells that relay chemosensory and 
mechanosensory feedback. When fitted to the experimen- 
tal data, the in silico model shows the existence of both 
stochastic and deterministic sources for variability in the 
respiratory rate. The upper and lower limits of variability 
seen experimentally are reproduced in silico using a range 
of network connection strengths. Simulations demonstrate 
that stochastic spiking in sensory relay neurons may be 
utilized to stabilize changes in variability when the respira- 
tory rate changes due to physiological demand. 
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